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We have investigated cosmological models with a self-interacting scalar field and a dissipative matter fluid 
as the sources of matter. Different variables are expressed in terms of a generating function. Exact solutions 
are obtained for one particular choice of the generating function The potential corresponding to this generating 
function is a standard tree-level potential arising in the perturbative regime in quantum field theory. With suitable 
choice of parameters, the scale factor in our model exhibits both inflationary behaviour in the early universe as 
well as an accelerating phase at late times with a decelerating period in between. It also satisfies the constraints 
for primeval nucleosynthesis and structure formation and seems to solve the cosmic coincidence problem. The 
solution exhibits a attractor nature towards a asymptotic de-sitter universe. 
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I. INTRODUCTION 

A number of recent observations ^ suggest that the Q m , the ratio of the matter density(baryonic+dark) to the critical density, 
is significantly less than unity suggesting that either the universe is open or that there is some other sources of this missing energy 
which makes tttotai ~ 1. The recent findings of BOOMERANG experiments [^] strongly suggests the second possibility of a 
flat universe. At the same time, the measurements of the luminosity -redshift relations observed for the 50 newly discovered type 
la supernova with redshift z > 0.35 |Q| indicate that at present the universe is expanding in an accelerated fashion suggesting a 
net negative pressure for the universe. 

Initial suggestions were to identify this missing energy density to a cosmological constant A [[jJBp. For a flat matter dominated 
universe with A in Einstein gravity observations strongly suggest ftj^ ~ 0.72. However, this possibility that A could be the 
dominant energy density has the drawback that the energy scale involved is lower than normal energy scale predicted by the 
most particle physics models by a factor of ~ 10~ 123 . An alternative source of energy density that may be admissible for this 
acceleration could be a dynamical A |Qj in the form of a scalar field with some self interacting potential If the energy density 
of this kind of source varies slowly with time, it mimics an effective cosmological constant. The idea of this candidate, called 
quintessence Q, is borrowed from the inflationary paradigm of the early universe. The difference, however, is that this new 
field evolves at a much lower energy scale. The energy density of this field, though dominant at present epoch, must remain 
sub-dominant at very early stages and should have evolved in such a way that it becomes comparable to the matter density p m 
now. When quintessence is modeled using a minimally coupled scalar field, in general, parameters need to be fine-tuned so as to 
ensure that p m and are of the same order today. This fine tuning problem has been termed as the cosmic coincidence problem. 
A new form of quintessence field called "tracker field' has been proposed to solve the cosmic coincidence problem. It has an 
equation of motion with an attractor like solution in the sense that for a wide range of initial conditions the equation of motion 
converges to the same solution. 

There are a number of quintessence models which have been suggested and most of these involve scalar fields with minimal 
coupling with potentials dominating over the kinetic energy of the field. A purely exponential potential is one of the widely 
studied cases [14]]. Inspite of the several advantages the energy density is not enough to make up for the missing part. In- 
verse power law is another form of the potential ( [^]- [[)]]) that has been considered extensively for quintessence models, in 
particular, for solving the cosmic coincidence problem. Though many of the problems are resolved successfully with this po- 
tential, the predicted value for the equation of state for the quintessence field, jq, is not in good agreement with the observed 
results. In search of suitable models that would eliminate the problems, new types of potentials, like Vb[cos hX<j) — l] p [ |lO| ] 
and Vq sin h(ay/koA^))^ [|5 1 1 ] have been considered, which have asymptotic forms like the inverse power law or exponential 



ones. Different physical considerations have lead to the study of other types of the potentials also |jj~2|]. Recently Saini et al jl3| ] 
have reconstructed the potential in context of general relativity and minimally coupled quintessence field from the expression of 
the luminosity distance cIl(z) as function of redshift obtained from the observational data. However, none of these potentials 
are entirely free of problems. Hence, there is still a need to identify appropriate potentials to explain current observations Jl4]]. 
Also it has been recently shown by Pietro and Demaret [jl5|] that for constant scalar field equation of state, which is a good 
approximation for a tracker field solutions, the field equations and the conservation equations strongly constrain the scalar field 
potential. Most of the widely used potential for quintessence, such as inverse power law one, exponential or the cosine form, are 
incompatible with these constraints. 

The CDM is in general considered to be a perfect fluid. However, in some scenarios, certain physical processes can make 
the CDM fluid effectively a dissipative one. In such a situation the fluid has an effective pressure that is negetive. Recently it 
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has been proposed that the CDM must self interact in order to explain the detailed structure of the galactic halos JT^]. This self 
interaction will create a viscous pressure whose magnitude will depend on the mean free path of the CDM particles. In a recent 
work Chimento et.al have shown that a mixture of minimally coupled self interacting scalar field and a perfect fluid is unable to 
drive the accelerated expansion as well as solve the cosmic coincidence problem at the same time JlTy. However, a mixture of a 
dissipative CDM with bulk viscosity and a minimally coupled self interacting scalar field can successfully achieve both features 
simultaneously. Also, as demonstrated in a recent paper by Zimdahl et. al. [|l8|] one can also have a negative n if there exists an 
interaction which does not conserve particle numbers. This may be due to the particle production out of gravitational field. In 
this case, the CDM is not a conventional dissipative fluid, but a perfect fluid with varying particle number. Substantial particle 
production is an event that occurs in the early universe. But Zimdahl et. al. have shown that even extremely small particle 
production rate can also cause the sufficiently negative ir to violate the strong energy condition. 

In this paper we have used a minimally coupled scalar field with a self interacting potential together with a matter fluid having 
a dissipative pressure over and above its positive equilibrium pressure. We have not assumed any particular model for this 
negative pressure. Instead, we have investigated what kind effects it has in the expansion of the universe. Unlike other works in 
scalar field cosmology with a dissipative pressure [ |l9| ] we have neither assumed the behaviour of the scale factor nor have we 
assumed any specific form of the potential. Rather we have expressed all the variables in terms of what we call the 'generating 



function' . For this we have followed the method described by Chimento et.al. [ 20 j with some additional assumptions. We 
have proceeded with a particular choice of the generating function for which the potential is constructed using a combination of 
different power-law functions of <f>. From the behaviour of decelerating parameter it has been shown that one can indeed generate 
both inflationary era in the early time and also an late time accelerating phase with a decelerating period in between. We have 
also investigated the stability and attractor structure of the general solutions of the field equations with this kind of potential and 
have found that for certain choices of the constants the solutions indeed exhibit attractor behaviour in the late times. 

II. FIELD EQUATIONS 

Let us consider a homogeneous, isotropic, spatially flat FRW universe with a line element 

3 

ds 2 = -dt 2 + a(t) 2 ^ dx h (!) 

i=l 

where a(t) is the scale factor. The energy density consists of a massive scalar field with a self-interacting potential V(cf>) 

S» u = 4>,^,v ~ SV [ \(t>, a (j)' a + V{4>) ) , (2) 
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together with a dissipative fluid having bulk viscosity as the only dissipative term: 

T^ v = pu^Uy + h^ v (p + 7r) , (3) 

where p is the energy density, p is the equilibrium pressure, tt is the bulk viscous pressure, and h^ v — g^ v 4- u^u u is the 
projection tensor. 

The independent field equations describing the system are 

4> + 3Hj>+^ = 0, (4) 

3H 2 = ^ 2 + V + p, (5) 
p + 3H(p + p + 7r) =0 (6) 

After some straightforward calculations one can construct the equation 

2ij = -0 2 - (p + p + TT). (7) 

For H 7^ 0, the system consists of three independent equations, viz., (^|), (^]) and (Q) and we have six unknowns, 
H, (f>,V(<f)), p,p and tt. Hence, three constraint equations are required in order to have a closed system of equations. We 
first assume that p,p and tt are related to <p by, 

p + p + n = m<p 2 , (8) 



2 



where m is an arbitrary constant. We emphasize that this choice is considered purely because it makes the system of equations 
simple to solve. We use this to see whether such a choice leads to any physically acceptable solutions. 

The nature of the potential determines the evolution of <j>. However, instead of choosing a particular form for V ((/)), we can 
alternatively describe the evolution of the scalar field by expressing as a function of cf> as, 



= n<p) 



(9) 



where F(<f) is called the generating function. In this approach, we start with a specific trajectory, <p — F (</>), in phase space 
for <fi and determine the potential that evolves the scalar field in that manner. However, choosing a potential V(<f>) or choosing a 
generating function, F(<j>) is not completely equivalent. Choosing a form for F((f>) is more restrictive than choosing a potential 
due to the following reason. A specific F(<f>) refers to a specific class of initial conditions. Hence the set of solutions represented 
by a specific form for F(4>) is a subset of the full set of solutions for the corresponding potential. Having constructed the 
potential function, we study the evolution of the system in phase space for a more general initial condition, i.e. ones which are 
not restricted by any specific form for F{(j>). 

Using equation (||), equation (^) and equation ([)]) one can write, 



dH 



= -(m + l)F(4>). 
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In this paper, we first choose F((f>) and then calculate the V{<p) that is consistent with our choice of F(<p) as well as with the 
Einstein's equations. The solutions for which the phas e spac e trajectory is given by = F(tfi) is a subset of the of the general 
class of solutions for this potential. (Later in subsection III C we will draw the phase space trajectory for a set of general solutions 
(i.e. with different initial conditions) including the one for which <p = F (</))) and check the conditions for this to be a stable 
solution.) 

For any given F((j>) which we term as the "generating function", one can, in principle, solve the system as follows: 



im = -^pl I F(<,Uln + tf„. 



p(<p) 

a(cf>) 



-3m 



ao exp 



H{cj>)F{cp)dcp 
F(4>) 



Po, 



V&)=3H 2 (cj>)--F(4>f 
p + 7T = mF((f))' 2 — p(<f>). 



p{4>), 



(ii) 

(12) 
(13) 

(14) 
(15) 



where Hq, po and ao are integration constant. With the assumption of an equation of state p = p(p), one can also calculate ir. 
Hence, our main aim is to choose F(4>) properly to have some physically acceptable behaviour for different variables. 



III. SOLUTIONS FOR F{<f>) = ujq 
A. Exact solutions 



With this choice of the generating function, the differential equation governing the time evolution of <j> is 

(J) = UHp. 

In this case, the exact solutions turn out to be, 



(j)(t) = cj) exp(wi), 



H{t) 



(m + 1) 



ui(f>Q exp(2o;i) + Hq 



ait) = ao exp 



$j(m + l) 



exp(2w£) + H Q t 



(16) 

(17) 
(18) 

(19) 



3 



p(t) 



(to + 1) 



u 4>o exp(2wt) 



-"0 \ 



(20) 



3(to + l)mui 2 <pQ cxp(4wi) 
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3TOwi?o0o exp(2wt) 



Po, 



f<m2 x , 3(to + 1V 2 4 c^ 2 



(21) 
(22) 



For an expanding Universe we must have H > 0. To ensure this, either rn < — 1 or uj < 0. We choose w < 0. It may be 
noted from equation (19) that the proper volume becomes zero at t — > — oo and hence this is taken to be the initial time for our 
model. The potential V(4>) given in equation (22) is the standard renormalizable tree level potential arising in the perturabative 
regime of quantum field theory. With u> < 0, the shape of the potential depends on the factor (u> + 3Hq). When Hq > — u>/3, it is 
the most simplified version of the potential for the hybrid inflation [pT|]. Similarly when Hq < — uj/3, it is the inverted quadratic 
potential and has been discussed by many authors for inflationary models (see pip and references therein). As has been pointed 
out before, the solution represented by equation ( jl7| ) is only one of the classes solutions (corresponding to the potential given in 
equation (p2]))for which the the initial values of <fi and </> are related by 4>i„itiai = ^"Pinitiai- 

Expressing the equation of state as, p = (7 — l)p, where 1 < 7 < 2, the bulk viscous pressure is given by, 



}(f>Q exp(2w£) ( u 



37771(771 + V 
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w 2 0q exp(4w7j) - 7p . 



(23) 



B. Behaviour of the solutions 

The behaviour of the deceleration parameter is shown in figure [l] The evolution of q(t) has the basic feature that we require 
for an acceptable form for the deceleration parameter. q(t) is negative to begin with resulting in an inflationary phase. It 
increases and subsequently q(t) becomes positive. At a later time, it drops below and finally saturates at a constant value 
below zero. This last feature of q(t) at late times when it becomes a negative constant produces the present day accelerated 
expansion. In fact it asymptotically becomes a de Sitter Universe. The behaviour of the equation of state for the scalar field 
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FIG. 1. Behaviour of the deceleration parameter with time 

w = (\4> 2 — V((f>))/(^4> 2 + V(4>)) in figure (Q) shows that at late times w = — 1 which essentially supports the existence of a 
cosmological constant in late time. 

An acceptable model should also satisfy constraints arising from cosmological nucleosynthesis and from structure formation. 
In order to keep these intact, the matter energy density should dominate over the energy density of the scalar field in the early 
universe. However, at late times the scalar field energy density should be more than that due to matter so that one can explain the 
missing energy of the universe. In figure (fy we have plotted the energy density for the matter (p) and for the scalar field (p^). 
We see the above constraints are satisfied. The matter energy density is dominant in the early time and so nucleosynthesis and 
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structure formation are unaffected. But at late time p^ decreases more slowly than p and hence ultimately it becomes greater 
than p This explains the missing energy associated with accelerated expansion of the universe. It can also be seen that although 
in early period the energy densities are of different order of magnitude but in late time they are of the same order. This gives 
a dynamical solution to the cosmic coincidence problem. Also at late times the ratio of these two energy densities becomes 
constant showing the tracking behaviour. We have also plotted the two density parameter Q m and f2^ in figure 4. 
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FIG. 2. Behaviour of w with time 




0.5 1 1.5 2 2.5 3 

t 



FIG. 3. Behaviour of p (dashed)and pj, (solid) with time 
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FIG. 4. Behaviour of Q m (dashed) and Q$ (solid) with time 
The variation of the viscous pressure ir is shown in g. 

C. Stability and attractor structure 

The solutions we discussed were those for which F(fa) — lj</>. In other words, for these set of solutions the initial conditions 
are constrained by the relation fa = ufa, where fa and fa are the initial values of <f> and <f>, respectively. These solutions are 
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FIG. 5. Behaviour of viscous pressure for 7 = 1 (dashed)and for 7 = 4/3 (solid), with time 

significant only if these are stable, i.e., only if more general initial conditions asymptotically evolve towards these phase space 
trajectories. 

This solution corresponds to the case when the expression for the Hubble parameter (H(t)), scale factor (<z(i)), deceleration 
parameter (q(t)) and the density (p(t)) given by equations (|l8|), (|l9|), ( ]20| ) and (pl"|). In this section we proceed to investigate 
whether or not the solutions described are stable and whether there are any constraints on the parameters for this. 

The general equation of motion for the potential given in equation ( p2[ ) is 

<i> = -3H4> - v V 3 + uj 2 (f> + 3H uj(t> (24) 
Together with the equation for the evolution of H{t), we can write the following set of three coupled differential equations, 

<j) = x 

x = -3Hx - 3 (ra + 1 ) CJ 2^3 + + 3HaLU( p 

J j = _(m±l) a;2 
2 

By equating the right hand sides of the equations to we obtain the following three critical points in the (</>, <j>) plane. 

The position of the critical points in the phase diagram changes with the value of H . Further, at oj = —3H all the three critical 
points merge. There is a transition in the nature of stability of the trajectories in (cj), 4>) space at this value of H . We discuss 
below the nature of the trajectories for different values of Hq. First of all, we should note that uj < 0. 

In figures (^) and ([7]), we have plotted the phase space trajectories for the scalar field with a variety of initial conditions that 
are not constrained by the particular choice of the generating function. Figure || corresponds to the case of Hq = 0.2. As we 
have chosen uj to be — 1, this is a case in which the potential has a local maximum. On the other hand, figure ^ corresponds to the 
case of Hq = 0.4 in which cf> = is a minimum for the potential. In the first case the phase space trajectories converge to either 
of the two non-zero values of <fi and is not an attractor like solution as the asymptotic nature of the solution depends on the initial 
conditions. But in the second case, however, the trajectories converge to (f> = 0, <j) — independent of the initial conditions and 
,the de-Sitter solution is an attractor for this case. 

rv. CONCLUSION 

The most important conclusion is that models with a self-interacting scalar field and a matter fluid having a negative pressure 
in addition to its positive equilibrium pressure can produce a scenario for the cosmological evolution in which one can have an 
inflationary phase to begin with, an accelerated phase at late times (like the present era) and a decelerating phase in-between. 
Recently Lopez and Matos [ p2[ | have shown that this kind of complete history for the scale factor can be described by a hyperbolic 
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FIG. 6. Phase space trajectories for Ho = 0.2 
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FIG. 7. Phase space trajectories for Ho = 0.4 

potential. But the physical origin of such potential is still not well known. But here we have shown that such kind of behaviour 
for the scale factor can be generated with a potential given in (22), which has been widely used by many authors for inflationary 
models. 

The behaviour of p and p^ in our model shows that although in early universe, p is greater than p^ which is necessary for 
different physical phenomena like nucleosynthesis and structure formation etc, in the late times, p^ starts dominating. This 
feature explains the missing energy density and also the ration of two energy densities becomes a constant in late time showing 
the "tracking nature". 

One should note that both the assumptions (8) and (16) play crucial role in our model. Given a barotropic equation of state 
between p and p one can not assume (8) and (16) at the same time if the dissipative pressure it is zero in our model as t hat will 
lead to an over-determined problem ( the number of unknowns will be less than the number of independent equations). Even if 
one assumes these two condition one can check that will lead to an negative equilibrium pressure which is not desirable. H ence 
the existence of dissipative pressure also plays an important role in our model. 

We have also studied the general equation of motion for the scalar field (equation (24)) for the potential (22) and have shown 
that for the choices of constant for which the potential is minimum at <\> = 0, the phase space diagram exhibit a attractor 
behaviour towards the asymptotic de-sitter solution. 

We want to mention that previously tracker and attractor solutions have been studied for scalar fields having inverse power 
law, exponential, cosine potential. But in all of these cases the equation of state it; is a constant in radiation era as well as in 
a matter dominated era. It was later shown by Pietro and Demaret [ |l5| ] that these kind of potential with a constant w is not 



7 



consistent with the field equations. In our case, the equation of state for the scalar field w is not a constant but it varies with the 
cosmic evolution and approaches towards -1 asymptotically showing the existence of a cosmological constant in late times. 
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